Emergence of zoonotic sporotrichosis in Brazil: a genomic epidemiology study.

BACKGROUND: Zoonotic sporotrichosis is a neglected fungal disease, whereby outbreaks are primarily driven by Sporothrix brasiliensis and linked to cat-to-human transmission. To understand the emergence and spread of sporotrichosis in Brazil, the epicentre of the current epidemic in South America, we aimed to conduct whole-genome sequencing (WGS) to describe the genomic epidemiology. METHODS: In this genomic epidemiology study, we included Sporothrix spp isolates from sporotrichosis cases from Brazil, Colombia, and the USA. We conducted WGS using Illumina NovaSeq on isolates collected by three laboratories in Brazil from humans and cats with sporotrichosis between 2013 and 2022. All isolates that were confirmed to be Sporothrix genus by internal transcribed spacer or beta-tubulin PCR sequencing were included in this study. We downloaded eight Sporothrix genome sequences from the National Center for Biotechnology Information (six from Brazil, two from Colombia). Three Sporothrix spp genome sequences from the USA were generated by the US Centers for Disease Control and Prevention as part of this study. We did phylogenetic analyses and correlated geographical and temporal case distribution with genotypic features of Sporothrix spp isolates. FINDINGS: 72 Sporothrix spp isolates from 55 human and 17 animal sporotrichosis cases were included: 67 (93%) were from Brazil, two (3%) from Colombia, and three (4%) from the USA. Cases spanned from 1999 to 2022. Most (61 [85%]) isolates were S brasiliensis, and all were reported from Brazil. Ten (14%) were Sporothrix schenckii and were reported from Brazil, USA, and Colombia. For S schenckii isolates, two distinct clades were observed wherein isolates clustered by geography. For S brasiliensis isolates, five clades separated by more than 100 000 single-nucleotide polymorphisms were observed. Among the five S brasiliensis clades, clades A and C contained isolates from both human and cat cases, and clade A contained isolates from six different states in Brazil. Compared with S brasiliensis isolates, larger genetic diversity was observed among S schenckii isolates from animal and human cases within a clade. INTERPRETATION: Our results suggest that the ongoing epidemic driven by S brasiliensis in Brazil represents several, independent emergence events followed by animal-to-animal and animal-to human transmission within and between Brazilian states. These results describe how S brasiliensis can emerge and spread within a country. FUNDING: Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, Brazil; the São Paulo Research Foundation; Productivity in Research fellowships by the National Council for Scientific and Technological Development, and Ministry of Science and Technology of Brazil.

Genomic epidemiology and antifungal-resistant characterization of Candida auris, Colombia, 2016-2021.

Since 2016, in Colombia, ongoing transmission of Candida auris has been reported in multiple cities. Here, we provide an updated description of C. auris genomic epidemiology and the dynamics of antifungal resistance in Colombia. We sequenced 99 isolates from C. auris cases with collection dates ranging from June 2016 to January 2021; the resulting sequences coupled with 103 previously generated sequences from C. auris cases were described in a phylogenetic analysis. All C. auris cases were clade IV. Of the 182 isolates with antifungal susceptibility data, 67 (37%) were resistant to fluconazole, and 39 (21%) were resistant to amphotericin B. Isolates predominately clustered by country except for 16 isolates from Bogotá, Colombia, which grouped with isolates from Venezuela. The largest cluster (N = 166 isolates) contained two subgroups. The first subgroup contained 26 isolates, mainly from César; of these, 85% (N = 22) were resistant to fluconazole. The second subgroup consisted of 47 isolates from the north coast; of these, 81% (N = 38) were resistant to amphotericin B. Mutations in the ERG11 and TAC1B genes were identified in fluconazole-resistant isolates. This work describes molecular mechanisms associated with C. auris antifungal resistance in Colombia. Overall, C. auris cases from different geographic locations in Colombia exhibited high genetic relatedness, suggesting continued transmission between cities since 2016. These findings also suggest a lack of or minimal introductions of different clades of C. auris into Colombia. IMPORTANCE: Candida auris is an emerging fungus that presents a serious global health threat and has caused multiple outbreaks in Colombia. This work discusses the likelihood of introductions and local transmission of C. auris and provides an updated description of the molecular mechanisms associated with antifungal resistance in Colombia. Efforts like this provide information about the evolving C. auris burden that could help guide public health strategies to control C. auris spread.

Bimodal distribution of azole susceptibility in Sporothrix brasiliensis isolates in Brazil.

Sporothrix brasiliensis is an emerging zoonotic fungal pathogen that can be difficult to treat. Antifungal susceptibility testing was performed on the mold phase of a convenience sample of 61 Sporothrix spp. isolates from human and cat sporotrichosis cases in Brazil using the Clinical and Laboratory Standards Institute standard M38. A bimodal distribution of azole susceptibility was observed with 50% (28/56) of S. brasiliensis isolates showing elevated itraconazole minimum inhibitory concentrations ≥16 µg/mL. Phylogenetic analysis found the in vitro resistant isolates were not clonal and were distributed across three different S. brasiliensis clades. Single nucleotide polymorphism (SNP) analysis was performed to identify potential mechanisms of in vitro resistance. Two of the 28 resistant isolates (MIC ≥16 mg/L) had a polymorphism in the cytochrome P450 gene, cyp51, corresponding to the well-known G448S substitution inducing azole resistance in Aspergillus fumigatus. SNPs corresponding to other known mechanisms of azole resistance were not identified in the remaining 26 in vitro resistant isolates.

Genomic description of human clinical Aspergillus fumigatus isolates, California, 2020.

Aspergillus fumigatus, an environmental mold, causes life-threatening infections. Studies on the phylogenetic structure of human clinical A. fumigatus isolates are limited. Here, we performed whole genome sequencing of 24 A. fumigatus isolates collected from 18 patients in U.S. healthcare facilities in California. Single-nucleotide polymorphism (SNP) differences between patient isolates ranged from 187 to 70 829 SNPs. For five patients with multiple isolates, we calculated the within-host diversities. Three patients had a within-host diversity that ranged from 4 to 10 SNPs and two patients ranged from 2 to 16 977 SNPs. Findings revealed highly diverse A. fumigatus strains among patients and two patterns of diversity for isolates that come from the same patient, low and extremely high diversity.

Aspergillus fumigatus is an environmental mold. It can cause a severe infection called aspergillosis in patients with weakened immune systems. We analyzed A. fumigatus DNA from patients at California hospitals. We described genetic diversity between samples from the same patients and among different patients. Our findings provide insight on using genomic sequencing to investigate aspergillosis in hospitals.

Investigation of a Candida auris outbreak in a skilled nursing facility - Virginia, United States, October 2020-June 2021.

Candida auris, an emerging multi-drug resistant organism, is an urgent public health threat. We report on a C. auris outbreak investigation at a Virginia ventilator skilled nursing facility. During October 2020-June 2021, we identified 28 cases among residents in the ventilator unit. Genomic evidence suggested ≥2 distinct C. auris introductions to the facility. We identified multiple infection and prevention control challenges, highlighting the importance of strengthening multi-drug resistant organism prevention efforts at ventilator skilled nursing facilities.

Improving the Annotation of the Venom Gland Transcriptome of Pamphobeteus verdolaga, Prospecting Novel Bioactive Peptides.

Spider venoms constitute a trove of novel peptides with biotechnological interest. Paucity of next-generation-sequencing (NGS) data generation has led to a description of less than 1% of these peptides. Increasing evidence supports the underestimation of the assembled genes a single transcriptome assembler can predict. Here, the transcriptome of the venom gland of the spider Pamphobeteus verdolaga was re-assembled, using three free access algorithms, Trinity, SOAPdenovo-Trans, and SPAdes, to obtain a more complete annotation. Assembler's performance was evaluated by contig number, N50, read representation on the assembly, and BUSCO's terms retrieval against the arthropod dataset. Out of all the assembled sequences with all software, 39.26% were common between the three assemblers, and 27.88% were uniquely assembled by Trinity, while 27.65% were uniquely assembled by SPAdes. The non-redundant merging of all three assemblies' output permitted the annotation of 9232 sequences, which was 23% more when compared to each software and 28% more when compared to the previous P. verdolaga annotation; moreover, the description of 65 novel theraphotoxins was possible. In the generation of data for non-model organisms, as well as in the search for novel peptides with biotechnological interest, it is highly recommended to employ at least two different transcriptome assemblers.

MycoSNP: A Portable Workflow for Performing Whole-Genome Sequencing Analysis of Candida auris.

Candida auris is an urgent public health threat characterized by high drug-resistant rates and rapid spread in healthcare settings worldwide. As part of the C. auris response, molecular surveillance has helped public health officials track the global spread and investigate local outbreaks. Here, we describe whole-genome sequencing analysis methods used for routine C. auris molecular surveillance in the United States; methods include reference selection, reference preparation, quality assessment and control of sequencing reads, read alignment, and single-nucleotide polymorphism calling and filtration. We also describe the newly developed pipeline MycoSNP, a portable workflow for performing whole-genome sequencing analysis of fungal organisms including C. auris.

Candida auris Whole-Genome Sequence Benchmark Dataset for Phylogenomic Pipelines.

Candida auris is a multidrug-resistant pathogen that represents a serious public health threat due to its rapid global emergence, increasing incidence of healthcare-associated outbreaks, and high rates of antifungal resistance. Whole-genome sequencing and genomic surveillance have the potential to bolster C. auris surveillance networks moving forward. Laboratories conducting genomic surveillance need to be able to compare analyses from various national and international surveillance partners to ensure that results are mutually trusted and understood. Therefore, we established an empirical outbreak benchmark dataset consisting of 23 C. auris genomes to help validate comparisons of genomic analyses and facilitate communication among surveillance networks. Our outbreak benchmark dataset represents a polyclonal phylogeny with three subclades. The genomes in this dataset are from well-vetted studies that are supported by multiple lines of evidence, which demonstrate that the whole-genome sequencing data, phylogenetic tree, and epidemiological data are all in agreement. This C. auris benchmark set allows for standardized comparisons of phylogenomic pipelines, ultimately promoting effective C. auris collaborations.

Hypertension and the roles of the 9p21.3 risk locus: Classic findings and new association data.

BACKGROUND: The band 9p21.3 contains an established genomic risk zone for cardiovascular disease (CVD). Since the initial 2007 Wellcome Trust Case Control Consortium study (WTCCC), the increased CVD risk associated with 9p21.3 has been confirmed by multiple studies in different continents. However, many years later there was still no confirmed report of a corresponding association of 9p21.3 with hypertension, a major CV risk factor, nor with blood pressure (BP). THEORY: In this contribution, we review the bipartite haplotype structure of the 9p21.3 risk locus: one block is devoid of protein-coding genes but contains the lead CVD risk SNPs, while the other block contains the first exon and regulatory DNA of the gene for the cell cycle inhibitor p15. We consider how findings from molecular biology offer possibilities of an involvement of p15 in hypertension etiology, with expression of the p15 gene modulated by genetic variation from within the 9p21.3 risk locus. RESULTS: We present original results from a Colombian study revealing moderate but persistent association signals for BP and hypertension within the classic 9p21.3 CVD risk locus. These SNPs are mostly confined to a 'hypertension island' that spans less than 60 kb and coincides with the p15 haplotype block. We find confirmation in data originating from much larger, recent European BP studies, albeit with opposite effect directions. CONCLUSION: Although more work will be needed to elucidate possible mechanisms, previous findings and new data prompt reconsidering the question of how variation in 9p21.3 might influence hypertension components of cardiovascular risk.

Mitochondrial Genome Sequences of the Emerging Fungal Pathogen Candida auris.

Candida auris is an emerging fungal pathogen capable of causing invasive infections in humans. Since its first appearance around 1996, it has been isolated in countries spanning five continents. C. auris is a yeast that has the potential to cause outbreaks in hospitals, can survive in adverse conditions, including dry surfaces and high temperatures, and has been frequently misidentified by traditional methods. Furthermore, strains have been identified that are resistant to two and even all three of the main classes of antifungals currently in use. Several nuclear genome assemblies of C. auris have been published representing different clades and continents, yet until recently, the mitochondrial genomes (mtDNA chromosomes) of this species and the closely related species of C. haemulonii, C. duobushaemulonii, and C. pseudohaemulonii had not been analyzed in depth. We used reads from PacBio and Illumina sequencing to obtain a de novo reference assembly of the mitochondrial genome of the C. auris clade I isolate B8441 from Pakistan. This assembly has a total size of 28.2 kb and contains 13 core protein-coding genes, 25 tRNAs and the 12S and 16S ribosomal subunits. We then performed a comparative analysis by aligning Illumina reads of 129 other isolates from South Asia, Japan, South Africa, and South America with the B8441 reference. The clades of the phylogenetic tree we obtained from the aligned mtDNA sequences were consistent with those derived from the nuclear genome. The mitochondrial genome revealed a generally low genetic variation within clades, although the South Asian clade displayed two sub-branches including strains from both Pakistan and India. In particular, the 86 isolates from Colombia and Venezuela had mtDNA sequences that were all identical at the base level, i.e., a single conserved haplotype or mitochondrial background that exhibited characteristic differences from the Pakistan reference isolate B8441, such as a unique 25-nt insert that may affect function.

Updates and Comparative Analysis of the Mitochondrial Genomes of Paracoccidioides spp. Using Oxford Nanopore MinION Sequencing.

The mitochondrial genome of the Paracoccidioides brasiliensis reference isolate Pb18 was first sequenced and described by Cardoso et al. (2007), as a circular genome with a size of 71.3 kb and containing 14 protein coding genes, 25 tRNAs, and the large and small subunits of ribosomal RNA. Later in 2011, Desjardins et al. (2011) obtained partial assemblies of mitochondrial genomes of P. lutzii (Pb01), P. americana (Pb03), and P. brasiliensis sensu stricto (Pb18), although with a size of only 43.1 kb for Pb18. Sequencing errors or other limitations resulting from earlier technologies, and the advantages of NGS (short and long reads), prompted us to improve and update the mtDNA sequences and annotations of two Paracoccidioides species. Using Oxford Nanopore and Illumina read sequencing, we generated high-quality complete de novo mitochondrial genome assemblies and annotations for P. brasiliensis (Pb18) and P. americana (Pb03). Both assemblies were characterized by an unusually long spacer or intron region (>50 kb) between exons 2 and 3 of the nad5 gene, which was moderately conserved between Pb03 and Pb18 but not similar to other reported sequences, except for an unassigned contig in the 2011 assembly of Pb03. The reliability of the insert missing from previous mtDNA genome assemblies was confirmed by inspection of the individual Nanopore read sequences containing nad5 coding DNA, and experimentally by PCR for Pb18. We propose that the insert may aid replication initiation and may be excised to produce a smaller structural variant. The updated mtDNA genomes should enable more accurate SNP and other comparative or evolutionary analyses and primer/probe designs. A comparative analysis of the mtDNA from 32 isolates of Paracoccidioides spp., using the SNPs of the aligned mitochondrial genomes, showed groupings within the brasiliensis species complex that were largely consistent with previous findings from only five mitochondrial loci.

The LUFS domain, its transcriptional regulator proteins, and drug resistance in the fungal pathogen Candida auris.

The LUFS domain (LUG/LUH, Flo8, single-strand DNA-binding protein [SSBP]) is a well-conserved and apparently ancient region found in diverse proteins and taxa. This domain, which has as its most obvious structural feature a series of three helices, has been identified in transcriptional regulator proteins of animals, plants, and fungi. Recently, in these pages (Wang et al., Protein Sci., 2019, 28:788-793), the first crystal structure of a LUFS domain was reported, for the human SSBP2, a transcriptional repressor. We briefly address how the new insights into LUFS structures might contribute to a better understanding of an important transcriptional activator of yeasts that contains the LUFS domain, Flo8, and consider how a focus on the LUFS domain and its variation could help us to understand etiologies of drug resistance in a recently emerged pathogenic fungus, Candida auris.

Molecular Epidemiology of Candida auris in Colombia Reveals a Highly Related, Countrywide Colonization With Regional Patterns in Amphotericin B Resistance.

Background: Candida auris is a multidrug-resistant yeast associated with hospital outbreaks worldwide. During 2015-2016, multiple outbreaks were reported in Colombia. We aimed to understand the extent of contamination in healthcare settings and to characterize the molecular epidemiology of C. auris in Colombia. Methods: We sampled patients, patient contacts, healthcare workers, and the environment in 4 hospitals with recent C. auris outbreaks. Using standardized protocols, people were swabbed at different body sites. Patient and procedure rooms were sectioned into 4 zones and surfaces were swabbed. We performed whole-genome sequencing (WGS) and antifungal susceptibility testing (AFST) on all isolates. Results: Seven of the 17 (41%) people swabbed were found to be colonized. Candida auris was isolated from 37 of 322 (11%) environmental samples. These were collected from a variety of items in all 4 zones. WGS and AFST revealed that although isolates were similar throughout the country, isolates from the northern region were genetically distinct and more resistant to amphotericin B (AmB) than the isolates from central Colombia. Four novel nonsynonymous mutations were found to be significantly associated with AmB resistance. Conclusions: Our results show that extensive C. auris contamination can occur and highlight the importance of adherence to appropriate infection control practices and disinfection strategies. Observed genetic diversity supports healthcare transmission and a recent expansion of C. auris within Colombia with divergent AmB susceptibility.

Draft Genome Sequences of Two Sporothrix schenckii Clinical Isolates Associated with Human Sporotrichosis in Colombia.

Sporothrix schenckii is a thermodimorphic fungal pathogen with a high genetic diversity. In this work, we present the assembly and similarity analysis of the whole-genome sequences of two clinical isolates from Colombia of S. schenckiisensu stricto.

Genome Diversity, Recombination, and Virulence across the Major Lineages of Paracoccidioides.

The Paracoccidioides genus includes two species of thermally dimorphic fungi that cause paracoccidioidomycosis, a neglected health-threatening human systemic mycosis endemic to Latin America. To examine the genome evolution and the diversity of Paracoccidioides spp., we conducted whole-genome sequencing of 31 isolates representing the phylogenetic, geographic, and ecological breadth of the genus. These samples included clinical, environmental and laboratory reference strains of the S1, PS2, PS3, and PS4 lineages of P. brasiliensis and also isolates of Paracoccidioides lutzii species. We completed the first annotated genome assemblies for the PS3 and PS4 lineages and found that gene order was highly conserved across the major lineages, with only a few chromosomal rearrangements. Comparing whole-genome assemblies of the major lineages with single-nucleotide polymorphisms (SNPs) predicted from the remaining 26 isolates, we identified a deep split of the S1 lineage into two clades we named S1a and S1b. We found evidence for greater genetic exchange between the S1b lineage and all other lineages; this may reflect the broad geographic range of S1b, which is often sympatric with the remaining, largely geographically isolated lineages. In addition, we found evidence of positive selection for the GP43 and PGA1 antigen genes and genes coding for other secreted proteins and proteases and lineage-specific loss-of-function mutations in cell wall and protease genes; these together may contribute to virulence and host immune response variation among natural isolates of Paracoccidioides spp. These insights into the recent evolutionary events highlight important differences between the lineages that could impact the distribution, pathogenicity, and ecology of Paracoccidioides. IMPORTANCE Characterization of genetic differences between lineages of the dimorphic human-pathogenic fungus Paracoccidioides can identify changes linked to important phenotypes and guide the development of new diagnostics and treatments. In this article, we compared genomes of 31 diverse isolates representing the major lineages of Paracoccidioides spp. and completed the first annotated genome sequences for the PS3 and PS4 lineages. We analyzed the population structure and characterized the genetic diversity among the lineages of Paracoccidioides, including a deep split of S1 into two lineages (S1a and S1b), and differentiated S1b, associated with most clinical cases, as the more highly recombining and diverse lineage. In addition, we found patterns of positive selection in surface proteins and secreted enzymes among the lineages, suggesting diversifying mechanisms of pathogenicity and adaptation across this species complex. These genetic differences suggest associations with the geographic range, pathogenicity, and ecological niches of Paracoccidioides lineages.

From NGS assembly challenges to instability of fungal mitochondrial genomes: A case study in genome complexity.

The presence of repetitive or non-unique DNA persisting over sizable regions of a eukaryotic genome can hinder the genome's successful de novo assembly from short reads: ambiguities in assigning genome locations to the non-unique subsequences can result in premature termination of contigs and thus overfragmented assemblies. Fungal mitochondrial (mtDNA) genomes are compact (typically less than 100 kb), yet often contain short non-unique sequences that can be shown to impede their successful de novo assembly in silico. Such repeats can also confuse processes in the cell in vivo. A well-studied example is ectopic (out-of-register, illegitimate) recombination associated with repeat pairs, which can lead to deletion of functionally important genes that are located between the repeats. Repeats that remain conserved over micro- or macroevolutionary timescales despite such risks may indicate functionally or structurally (e.g., for replication) important regions. This principle could form the basis of a mining strategy for accelerating discovery of function in genome sequences. We present here our screening of a sample of 11 fully sequenced fungal mitochondrial genomes by observing where exact k-mer repeats occurred several times; initial analyses motivated us to focus on 17-mers occurring more than three times. Based on the diverse repeats we observe, we propose that such screening may serve as an efficient expedient for gaining a rapid but representative first insight into the repeat landscapes of sparsely characterized mitochondrial chromosomes. Our matching of the flagged repeats to previously reported regions of interest supports the idea that systems of persisting, non-trivial repeats in genomes can often highlight features meriting further attention.

The Dynamic Genome and Transcriptome of the Human Fungal Pathogen Blastomyces and Close Relative Emmonsia.

Three closely related thermally dimorphic pathogens are causal agents of major fungal diseases affecting humans in the Americas: blastomycosis, histoplasmosis and paracoccidioidomycosis. Here we report the genome sequence and analysis of four strains of the etiological agent of blastomycosis, Blastomyces, and two species of the related genus Emmonsia, typically pathogens of small mammals. Compared to related species, Blastomyces genomes are highly expanded, with long, often sharply demarcated tracts of low GC-content sequence. These GC-poor isochore-like regions are enriched for gypsy elements, are variable in total size between isolates, and are least expanded in the avirulent B. dermatitidis strain ER-3 as compared with the virulent B. gilchristii strain SLH14081. The lack of similar regions in related species suggests these isochore-like regions originated recently in the ancestor of the Blastomyces lineage. While gene content is highly conserved between Blastomyces and related fungi, we identified changes in copy number of genes potentially involved in host interaction, including proteases and characterized antigens. In addition, we studied gene expression changes of B. dermatitidis during the interaction of the infectious yeast form with macrophages and in a mouse model. Both experiments highlight a strong antioxidant defense response in Blastomyces, and upregulation of dioxygenases in vivo suggests that dioxide produced by antioxidants may be further utilized for amino acid metabolism. We identify a number of functional categories upregulated exclusively in vivo, such as secreted proteins, zinc acquisition proteins, and cysteine and tryptophan metabolism, which may include critical virulence factors missed before in in vitro studies. Across the dimorphic fungi, loss of certain zinc acquisition genes and differences in amino acid metabolism suggest unique adaptations of Blastomyces to its host environment. These results reveal the dynamics of genome evolution and of factors contributing to virulence in Blastomyces.

Genome update of the dimorphic human pathogenic fungi causing paracoccidioidomycosis.

Paracoccidiodomycosis (PCM) is a clinically important fungal disease that can acquire serious systemic forms and is caused by the thermodimorphic fungal Paracoccidioides spp. PCM is a tropical disease that is endemic in Latin America, where up to ten million people are infected; 80% of reported cases occur in Brazil, followed by Colombia and Venezuela. To enable genomic studies and to better characterize the pathogenesis of this dimorphic fungus, two reference strains of P. brasiliensis (Pb03, Pb18) and one strain of P. lutzii (Pb01) were sequenced [1]. While the initial draft assemblies were accurate in large scale structure and had high overall base quality, the sequences had frequent small scale defects such as poor quality stretches, unknown bases (N's), and artifactual deletions or nucleotide duplications, all of which caused larger scale errors in predicted gene structures. Since assembly consensus errors can now be addressed using next generation sequencing (NGS) in combination with recent methods allowing systematic assembly improvement, we re-sequenced the three reference strains of Paracoccidioides spp. using Illumina technology. We utilized the high sequencing depth to re-evaluate and improve the original assemblies generated from Sanger sequence reads, and obtained more complete and accurate reference assemblies. The new assemblies led to improved transcript predictions for the vast majority of genes of these reference strains, and often substantially corrected gene structures. These include several genes that are central to virulence or expressed during the pathogenic yeast stage in Paracoccidioides and other fungi, such as HSP90, RYP1-3, BAD1, catalase B, alpha-1,3-glucan synthase and the beta glucan synthase target gene FKS1. The improvement and validation of these reference sequences will now allow more accurate genome-based analyses. To our knowledge, this is one of the first reports of a fully automated and quality-assessed upgrade of a genome assembly and annotation for a non-model fungus.

The complex task of choosing a de novo assembly: lessons from fungal genomes.

Selecting the values of parameters used by de novo genomic assembly programs, or choosing an optimal de novo assembly from several runs obtained with different parameters or programs, are tasks that can require complex decision-making. A key parameter that must be supplied to typical next generation sequencing (NGS) assemblers is the k-mer length, i.e., the word size that determines which de Bruijn graph the program should map out and use. The topic of assembly selection criteria was recently revisited in the Assemblathon 2 study (Bradnam et al., 2013). Although no clear message was delivered with regard to optimal k-mer lengths, it was shown with examples that it is sometimes important to decide if one is most interested in optimizing the sequences of protein-coding genes (the gene space) or in optimizing the whole genome sequence including the intergenic DNA, as what is best for one criterion may not be best for the other. In the present study, our aim was to better understand how the assembly of unicellular fungi (which are typically intermediate in size and complexity between prokaryotes and metazoan eukaryotes) can change as one varies the k-mer values over a wide range. We used two different de novo assembly programs (SOAPdenovo2 and ABySS), and simple assembly metrics that also focused on success in assembling the gene space and repetitive elements. A recent increase in Illumina read length to around 150 bp allowed us to attempt de novo assemblies with a larger range of k-mers, up to 127 bp. We applied these methods to Illumina paired-end sequencing read sets of fungal strains of Paracoccidioides brasiliensis and other species. By visualizing the results in simple plots, we were able to track the effect of changing k-mer size and assembly program, and to demonstrate how such plots can readily reveal discontinuities or other unexpected characteristics that assembly programs can present in practice, especially when they are used in a traditional molecular microbiology laboratory with a 'genomics corner'. Here we propose and apply a component of a first pass validation methodology for benchmarking and understanding fungal genome de novo assembly processes.